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Abstract 

While  it  is  common  knowledge  that  ion  beams  are  easily  neutralized  for  both  current  and  charge  density  using  a 
variety  of  means,  the  precise  process  of  neutralization  remains  unknown.  With  the  increasing  importance  of  electric 
propulsion,  and  in  particular  micropropulsion  systems,  this  question  is  of  significant  importance.  Additionally,  it 
has  bearing  on  thruster  design,  space  instrument  calibration,  electrodynamic  tethers,  and  ionospheric  research.  A 
review  of  the  present  state  of  knowledge  on  this  topic  is  presented  as  well  as  results  from  ion  beam  simulations 
using  2D  and  3D  Particle-in-Cell  (PIC)  codes.  We  investigate  both  the  early  ‘Tilling”  problem  of  the  beam  starting 
to  move  away  from  the  spacecraft  and  the  steady  state  problem  where  the  beam  encounters  a  wall  at  an  infinite 
distance  from  the  spacecraft. 

PACS  codes: 

52.65.Rr  Particle-in-Cell  Method 
52.75.Di  Accelerators  and  Propulsion 
52.40.-w  Plasma  interactions 
52.25.-b  Plasma  Properties 


1.  Introduction 

Ion  beam  neutralization  during  operation  of  electric  propulsion  devices  requires  both  current  and  charge  density 
matching  with  an  emitted  electron  beam.  We  describe  the  phenomenon  as  “current  coupling”  and  note  that  it  is 
easily  accomplished  in  practice,  yet  the  exact  process  remains  unknown.  Currently  the  neutralization  process  is 
described  through  an  “effective  collision  frequency”  that  binds  electrons  to  the  ion  beam.  As  electric  propulsion 
becomes  more  prevalent  in  space  missions,  this  question  is  of  significant  importance.  Proper  modeling  of  current 
coupling  and  neutralization  will  enable  development  of  low-current  neutralizers  and  optimization  of  neutralizers  for 
micropropulsion  devices.  Explanation  of  the  effective  collision  frequency  also  has  bearing  on  space  instrument 
calibration,  electrodynamic  tethers,  and  ionospheric  research. 

In  the  early  years  of  electric  propulsion  research,  the  ion  beam  neutralization  question  was  one  of  the 
fundamental  issues  for  successful  development  of  this  promising  technology.  A  dense  ion  beam  requires  space 
charge  neutralization  to  avoid  a  potential  barrier  that  can  divert  or  reflect  the  beam.  The  vehicle  on  which  the 
thruster  operates  needs  current-neutrality  to  avoid  excessive  charging.  In  the  context  of  collisionless  plasma  theory, 
achieving  both  current  and  charge  neutrality  with  the  same  source  of  electrons  appears  to  be  nearly  impossible 
owing  mostly  to  the  large  difference  in  mass  between  electrons  and  the  ions.  For  example,  define  the  ion  flux, 
Fi  =  N ivi  and  the  net  electron  flux,  Fe  =  }iNevet  where  N  is  density,  v.  is  the  ion  drift  velocity,  and  vet  is  the 
electron  thermal  velocity  for  a  simple  effusion  model.  Achieving  both  equal  density  and  flux  requires  vet  =  Avi .  For 
example,  a  1-keV  Xenon  beam  has  ^=38,000  m/s  so  a  matching  electron  velocity  requires  a  source  temperature  of 
about  0.05  eV;  a  challenging,  but  not  impossible  number.  Replacing  this  simple  electron  effusion  assumption  with 
an  injected  distribution  does  not  change  the  misleading  picture  of  a  precise  balance  requirement.  A  higher 
temperature,  lower  density  electron  source  will  lead  to  a  positive  potential  well  that  does  trap  electrons,  but  then  the 
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theory  must  explain  by  what  process  the  trapped  electrons  shed  energy  so  as  to  actually  fill  the  well  Another 
observation  is  that  when  ion  beams  and  neutralizers  are  operated  in  conducting  vacuum  tanks  the  currents  are 
closely  coupled  even  though  the  grounding  tank  eliminates  the  charge  accumulation  that  could  provide  feedback  for 
current  balance  so  it  appears  that  one  or  more  plasma  mechanisms  must  be  responsible  for  curreS  couplbg 

Possibly  first  pointed  out  by  L.  Spitzer  in  1952  [uncited  note  in  Seitz  et  al.  1961],  electric  propulsion  plumes 

tests' th  bC  PI0Per  y  electrons  or  else  severe  space-charge  effects  would  result.  Before  the  first  space 

?r; ‘h,  .e  fn°r  d0UbtS  .to  the  stablhty  of  any  neutralization  approach  to  the  ion  beam  created  by  an 
i  ^  fenT  ldea was  for  ^tralization  to  occur  shortly  after  emission  to  prevent  beam  refom 
However,  there  was  lack  of  understanding  as  to  how  the  electrons  would  stay  within  the  beam  if  they  were  injected 
”  l  the  neutralization  process  was  unstable  to  small  perturbations.  Failing  to  properly  neutralize  the  beam  would 
cause  a  dramatic  reduction  in  thrust,  as  a  significant  portion  of  the  beam  would  return  to  the  spacecraft  This 
problem  was  first  addressed  by  the  Ramo-Wooldridge  staff  in  their  review  of  electrostatic  propulsion  in  1960  1 

““ enough  to  provide  a  satisfact°^  -  -  * 

uT  *5?  ueXt  fT  rarS’many  theorists  who  looked  at  U  and  2-D  models  predicted  growing  instabilities  that 

Freii^findinz  Sat'osedi  T  S°me  findblgs  Pointed  towards  the  Possibility  of  neutralization,  such  as 

,  h  .  finding  that  oscillations  m  electron  current  density  aided  in  neutralization  or  Mirels,  (1961)3  who  found 

that  emitter  location  and  thermal  electrons  motion  would  not  significantly  impact  beam  neutralization  Others 
pointed  towards  potential  problems,  such  as  Seitz  et  al.4  who  noted  that,  if  electrons  were  left  to  drift  into  the 
exhaust  beam  from  a  stationary  start  outside  of  it,  the  center  of  the  beam  would  still  develop  a  virtual  cathode  and 
reduction.  Some  of  the  earliest  computational  studies  were  brought  to  bear  on  the  problem,  and 
and  Kooyers  >  asmf  a  one-dimensional  PIC  code  in  1963  were  able  to  provide  a  neutralized  beam  when 
electrons  were  injected  at  velocities  lower  than  the  directed  ion  velocity.  Fluctuations  in  the  space  charge  field 

showing  ^  b^m'  ,TW°  years  later  Wadhwa  et  al.6  performed  a  two-dimensional  PIC  simulation 

nntThS  1 6  eC,tr°nS  WOa  d  oacllJate  Wltlun  tbe  beam  to  allow  for  neutralization,  but  theorized  that  oscillations  was 
not  the  only  mechanism  at  work.  One  method  suggested  was  that  fluctuations  in  the  space-charge  field  allowed  for 
en  opy  increase  to  mix  the  electrons,  but  these  fluctuations  were  not  found  downstream  of  the  neutralizer 

from  *  Spa“EleCtriC  R°cket  Test  I  (SERT  I)  found  that  it  was  quite  easy  to  neutralize  ion  beams  in  space 

from  a  simple  neutralizer  geometry.  In  a  senes  of  tests  it  was  shown  that  the  ion  thruster  developed  thrust  at  a  level 
indicating  complete  beam  neutralization.  After  SERT  I,  proof  of  concept  was  achieved  and  the  theoretical 
discussion  of  beam  neutralization  was  dropped  in  favor  of  engineering  new  thrusters.  Studies  after  SERT  I  include 

fJaUat.10nS  °foneAUtya  ,Zer  P  aCCTnt  °Ptimizatlon  of  the  thrusters,  and  simulations  to  analyze  spacecraft-plume 
J n ™ M  a  A  ,  7  “ncal  simulations  of  neutralization  have  been  performed  recently,  including  Othmer  et 

?!' Usmg  a  relativistic  3D  PIC  simulation  and  Tajmar  and  Wang  investigating  FEEP  neutralizer  placement.12 
Neither  proposed  a  solution  to  the  mechanism  of  current  coupling. 

wwS*6  deC,afS  °f  ieSearCb  3nd  the  imPlementation  of  electric  propulsion  devices,  the  detailed  process  by 
which  an  ionized  beam  is  neutralized  m  space  is  still  unknown.  Assorted  methods  to  fit  data  with  theory  have  been 
found,  but  the  actual  process  has  yet  to  be  studied  in  sufficient  detail  to  folly  understand  the  subject.  Further  new 
electric  micropropulsion  devices  such  as  the  FEEP  or  the  colloidal  thrusters  or  large  arrays  of  ion  and  Hall  thmsters 
are  still  not  guaranteed  to  behave.  We  might  also  desire  a  means  to  predict  and  optimize  neutralizer  operations. 
Thus,  a  simulation  technique  exhibiting  beam  coupling  is  needed.  Additionally,  results  from  ion  beam 
neutralization  modeling  will  be  applicable  to  ion  beams  for  instrument  calibration,  electrodynamic  tethers 
ionospheric  research,  and  fundamental  plasma  physics.  J 

tr.^Urtimmedia?t8°aliS  t0  dfenTe  if  a  standard  Particle-In-Cell,  PIC,  technique  is  adequate  or  if  additional 
understand  a«d  capture  the  current  coupling  process.  We  present  a  series  of  simulations  using 

,  r  f  d  K  3S  Wel  aS  f  30  P  C/DSMC  code-  ’  These  simulations  show  the  dependence  of  the  beam 
neutralization  on  beam  energy  and  neutralization  current  during  the  initial  “filling”  problem  and  the  lack  of  coupling 

™du”s'„af  £  3  D1 "S  Paper  alS°  “  “  °f  Valida,i0n  0f  ’ PIC- 


2.  PIC  Simulations  and  Discussion 


2-D  simulations  centered  first  on  examining  the  dependence  of  the  problem  on  assorted  parameters  such  as 
current  density,  injection  energy,  and  thermal  energy.  They  also  gave  a  feel  for  how  the  system  “should”  respond  in 
a  normal  PIC  simulation.  We  examined  current-matched  and  density  matched  ion  and  electron  beams  at  different 


* 


energies.  The  domain  used  for  2-D  simulation  is  shown  in  Figure  1.  In  the  current-matched  cases,  the  beam  stalled 
while  in  the  density-matched  cases  the  beam  propagated  without  hindrance.  Figures  2  and  3  show  typical  vx-x  phase 
space  plots  for  these  simulations.  Parameters  for  the  current-matched  and  density-matched  cases  are  shown  in  Table 
1.  Additional  study  was  performed  examining  the  ion/electron  speed  ratio.  This  highlighted  that  while  faster 
electrons  were  less  likely  to  neutralize  the  beam  by  remaining  in  the  potential  well  created  by  the  ions,  it  was  not  the 
primary  factor.  Higher  energy  beams  (1000  eV)  neutralized  readily  even  with  electrons  moving  much  faster  while 
low  energy  beams  (100  eV)  could  not  be  brought  to  neutralize  without  an  excess  of  electrons  (density  matching). 
Results  can  be  seen  in  Table  2. 

In  all  cases  the  majority  of  electrons  were  trapped  in  the  potential  well  created  by  the  ion  beam.  In  current- 
matched  cases,  the  electrons  were  unable  to  fully  neutralize  the  beam  and  were  accelerated  to  higher  energies  by  the 
resulting  space-charge.  The  electrons  oscillated  along  the  length  of  the  beam,  and  presumably  along  other 
dimensions  as  well,  creating  a  fairly  uniform  distribution  along  the  beam.  In  cases  of  beam  stalling,  a  distinct  “inner 
core”  of  electrons  can  clearly  be  seen  on  the  phase  space  plots,  as  in  Figure  2.  The  front  edge  was  located  at  the 
front  of  the  ion  beam  where  the  particles  were  accelerated  due  to  a  “surfing”  effect  from  the  potential  structure.  The 
trailing  edge  was  located  along  the  point  of  highest  potential  in  the  beam,  indicating  that  electrons  did  not  fully 
return  to  the  injection  plane  to  neutralize  the  space  charge  that  develops  there. 

However,  one  can  argue  that  this  filling  problem  is  not  what  the  thruster  sees  at  all,  as  eventually  the  end  of  the 
beam  will  fully  couple  with  the  background  plasma  and  decouple  from  the  spacecraft.  It  is  certainly  not  what  is  seen 
after  more  that  a  few  milliseconds  in  a  vacuum  chamber  as  the  beam  will  impinge  upon  the  wall  of  the  chamber  or 
whatever  baffling  mechanisms  are  in  place.  Therefore,  the  beam  would  have  definite  ends  on  both  sides.  This 
makes  the  bounded  problem  capable  of  analyzing  which  side  the  electrons  exit  the  problem:  we  can  determine  if 
there  is  any  sort  of  ion  drag  that  indicates  current  coupling. 

To  investigate  current  coupling  six  simulations  were  performed,  varying  between  frozen  and  unfrozen  ions  and 
electron  currents  at  0.9,  1 .0,  and  1 . 1  times  the  ion  current.  The  base  was  the  unfrozen  case  at  90%  of  the  ion  current. 
Spreading  experienced  by  the  unfrozen  ions  was  the  sole  observable  difference  between  the  simulations.  The 
variation  in  electron  current  had  no  significant  impact  on  the  potential  of  the  beam,  and  all  maintained  a  potential  of 
approximately  3-5%  of  beam  energy. 

Using  the  base  case,  the  total  current  collected  at  each  end  of  the  simulation  was  monitored  by  the  code.  If  the 
electrons  were  coupling  with  the  ions,  there  would  be  a  distinct  difference  in  the  current  collected  at  each  end  as 
well  as  a  difference  between  the  current  collected  when  the  ions  were  frozen  as  opposed  to  unfrozen.  Four 
simulations  investigating  this  scenario  were  performed,  as  outlined  in  Table  4.  We  initially  found  that  there  was  a 
predisposition  towards  electrons  leaving  through  the  right  hand  side,  but  this  may  have  been  caused  by  ballistic 
electrons  that  begin  with  enough  energy  to  move  across  the  beam  and  out  the  other  side  without  slowing  and 
reflecting.  To  eliminate  this,  we  began  injecting  electrons  on  both  sides,  each  with  half  the  total  electron  current. 

To  determine  the  fraction  of  electrons  leaving  the  domain  on  each  side,  one  can  develop  a  simple  formula. 
From  current  conservation,  the  current  passing  through,  say,  the  right  electrode  is  equal  to  the  current  from  the 
electron  and  ion  sources  as  such: 

h  =  +  +  (1) 

where  the  subscripts  denote  the  ion  or  electron  current  with  i  and  e  respectively  and  the  collection  side  for  R. 
Superscripts  denote  the  injection  side.  Assuming  that  the  fraction  of  electrons /leaving  each  side  is  equivalent  for 
particles  emitted  from  either  side,  we  can  simplify  (1)  to 

4=-wf+/«/r  (2) 

where  If  is  the  total  electron  current  emitted  from  both  sides.  Since  if  is  a  component  of  If ,  we  can  define 
it  as  If  =  (3 If ,  so  p  is  the  fraction  of  electrons  emitted  from  the  right  side.  Similarly,  the  electron  current  is  some 
fraction  a  of  the  ion  current,  so  we  can  define  If  =  al{ .  This  allows  us  to  rewrite  (2)  as 


Solving  for  fr ,  we  find 


IR=-It-aIt(l3-£). 


(3) 


Similarly,  we  find 


Oil. 


fi  —  -  J  +  7 
aL 


(4) 


(5) 


where  fR  +  fL-  1  and  /?  +  7  =  1 .  A  similar  solution  for  fR  and  fL  can  be  found  in  the  case  of  frozen  ions  where 
no  ion  current  is  collected.  Full  current  coupling  would  be  indicated  by  fR  =  1  as  all  the  electrons  follow  the  ions 

ou  ofthe  problem  on  the  right  side.  If  there  is  no  preference  for  either  side,  the  ratio  fR/fL  should  scale  as  the  ratio 

CHUeCtmg  C,UITent'  ,The  current  on  the  electrode  was  calculated  and  recorded  over  the  duration  ofthe 

!imu“0n  Pai0i’  a  brief  Wind0W  *•  before  to 

In  XPDP2,  the  current  diagnostic  is  related  to  the  electrode  at  each  end,  which  also  injects  particles  into  the 
th  “  This  requirement  forced  us  to  perform  simulations  in  both  with  equal  area  electrodes  at  each  end  and  where 

beam8  Thk'llirt  5  the  S1f  °f  the  left  hand  electrode  to  accommodate  the  spreading  ofthe  ion 

beam  This  affects  the  collected  currents,  as  shown  above,  but  does  not  invalidate  the  simulations 

H  -  OsToa"  44  L  >W~  f’T?  n°  f»  the  diection  of  electron  leaving  the  simulation  when 

aS/«-/‘-(l5-  Ttere  s,”'tls  “  be  a  slight  preference  for  drift  out  the  right  in  case  all  cases 

SLSTuI  tof  5f.preferf‘“  K  **  i"  electrode  sizes  as  it  is  not  affected  by  freezing  or 

Tkrlt  7  conclusively  demonstrates  the  lack  of  cinrent  coupling  in  steady-state  PIC  simulations 
The  role  of  perveance  m  beam  neutralization  remains  under  investigation.  As  perveance  decreased  in  the 
beams  the  ability  to  neutralize  the  beam  dramatically  increased.  Also  looking  *  the  Tam^anance  in 
susceptibility  to  electron  velocity,  perveance  likely  plays  a  large  role  in  neutralization  stability. 

, lhe  3-D  domam  shown  m  Figure  4  was  generated  similar  to  2-D  and  the  filling  problem  simulations  were 

took’Ti  USC  °f  \3'D  C,°de  t0  m°del  this  essentia,1y  2-D  problem  was  performed  as  a  trial  before  simulations 
on  the  more  physical  case  with  a  separate,  localized  neutralizer.  The  results  obtained  match  the  2  D 
simulations  closely.  A  potential  well  develops  as  in  the  2-D  case,  as  seen  in  Figure  5.  The  ion  phases  becomes 
more  rounded  as  alternate  routes  around  the  potential  wells  are  found  as  in  Figures  6  and  7.  Ions  accelerated  out 
o  e  well  are  also  seen  in  3-D.  A  noticeable  effect  of  downstream  Neumann  (floating)  boundary  conditions  was 
lscanSeCaoid  e  y  mje,CIf  electrons>  which  could  easily  reach  the  downstream  boundary  and  were  allowed  to 

effectas  Senii^D1 8  ThscorTn6  mT  Y  deePenifg  Potential  well  thus  creating  a  vYry  visible  “inner  core” 

1  C°re  C0U,ld be  seen  t0  all8n  lts  leading  edge  with  the  accelerated  ion  front  and  the  trailing 

edge  with  the  pomt  of  beam  reversal,  again  as  seen  in  2-D.  trailing 

3.  Conclusions 

The  inability  of  the  steady  state  problem  to  demonstrate  current  coupling  indicates  that  there  are  mechanisms 
tappemng  to,  a,,  no,  considered  in  standard  PIC  codes.  While  i,  can  be  shown  to,  to  c™„  l«h” 
order  to  provide  a  stable  start  to  the  ion  beam  during  the  initial  “filling”  problem,  by  the  time  the  problem  reaches 

that  the  foil  nhvsi«n^ChamST  ***  deCtr°nS  down  the  ion  beam  as  simulated  in  standard  PIC.  It  is  clear 
nmJht  ? UFu  y  ♦  °fbeatn  neutrallzatl°n  are  not  understood  at  present.  Previous  work  on  the  subject  has  been 

undpr?dH  C  earuanSWer'  FVrther  devel°Pment  of  models  and  continued  study  is  necessary  for  a  complete 

understanding  of  the  ion  beam  neutralization  problem.  y  P 

This  series  of  2D  simulations  has  demonstrated  that  what  can  be  called  a  standard  approach  to  explicit  PIC  does 
demonstrate  charge  neutralization,  but  not  current  neutralization.  Since  the  charge  neutralization  was  present  in  all 
die  runs  including  those  with  frozen  ions  we  may  conclude  that  the  electrons  “fall”  into  the  positive  potential  well  by 
exchanging  energy  with  other  electrons  via  the  action  ofthe  electrostatic  field.  Since  the  simulations  do  not  exhibh 

ZZZf ft/*  may  alS0  COnClUfe  &at  m0ti0n  0f  the  ionS  is  n0t  any  coherent  phase  velocity  So  tte 

ostatic  field  i.e.  propagating  electrostatic  waves  which  would  entrain  electrons.  Since  these  waves  are 

weT  ri8  r? the  “?  bef,m  at  a  Veloclty  sma11  comPared  to  the  electron  velocity,  the  resonant  coupling  should  be 
ak.  Our  future  work  will  purposely  mject  ion  density  waves  to  determine  what  is  required  for  PIC  to  exhibit  this 

“  “7  T »PPtos«,  numerically.  If,  however,  we  cm,  «2s,  tose  sfmtodotte  fan 
formulate  the  following  hypothesis:  either  an  ion  beam  couples  to  its  neutralizer  electrons  due  to  correlations  or 
structures  in  the  ion  population  that  can  mduce  a  down  stream  phase  velocity  in  the  electrostatic  field  an  effect  that 

Sfoat  eleCtr°nS  are  reSPOnding  t0  SOmething  “»  Coul0mb  COlta  -  a  sub-Debye 
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Constant  Simulation  Parameters  1 

ncx 

ncy 

4 

4 

r 

elec 

128 

128 

lifiH 

1.0E-11 

0.025 

n, 

_ * _ i 

W 

mion 

wm 

4 

1E11 

le4 

2.18e-25 

0 

0.1 

0.01 

Simulation  Variables  and  Results 

Case 

3i 

i/ji 

E  . 

XI 

^raax 

P 

i.i 

1.21E+19 

1.00 

100 

111.8380 

4.85E-7 

1.2 

1.21E+19 

3.46 

100 

1.7708 

4.85E-7 

1.3 

3.83E+19 

1.00 

1000 

20.6984 

4.85E-8 

1.4 

3.83E+19 

1.09 

1000 

1.5144 

4.85E-8 

Table  1:  2-D  Simulation  parameters  and  current-density  comparison  results 


Case 

Base 

case 

Te  (eV) 

ve 

(m/s) 

vx 

Virtual 

Cathode? 

2.1 

1.5 

1.0 

419386 

10.94 

Yes 

2.2 

1.5 

0.1 

132621 

3.46 

No 

2.3 

1.3 

8.36E-4 

12123 

1.0 

Yes 

2.4 

1.3 

2.09E-4 

6061.5 

0.5 

Yes 

Table  2:  2-D  velocity  ratio  simulation  parameters  and  results 


Case 

Ions 

p 

a 

4L  [V] 

3.1 

Frozen 

0 

0.9 

0.38 

3.2 

Warm 

0 

0.9 

0.59 

3.3 

Frozen 

0 

1.0 

0.69 

3.4 

Warm 

0 

1.0 

1.73 

3.5 

Frozen 

0 

1.1 

0.50 

3.6 

Warm 

0 

1.1 

1.53 

Table  3:  Steady  state  simulation  parameters  and  results.  “Warm”  ions  indicate  the  ions  were  mobile  and  allowed  to 
develop  from  a  cold  beam.  “Frozen”  ions  indicate  particles  were  place  in  the  problem  by  the  code  in  a  perfectly 
cold  beam. 


Simulation  # 

Ions 

p 

a 

R.  Electrode  Size 

fl 

/fl 

4.1 

Warm 

A 

0.9 

2 

0.33 

0.59 

4.2 

Cold 

! A 

0.9 

2 

0.35 

0.61 

4.3 

Warm 

XA 

1.0 

2 

0.32 

0.67 

4.4 

Frozen 

A 

1.0 

1 

0.46 

0.49 

Table  4:  Current  coupling  simulation  parameters  and  results.  “Warm”  ions  indicate  the  ions  were  mobile  and 
started  from  the  developed  case  3.2.  “Cold”  ions  indicate  the  ions  were  immobile,  but  started  from  the  developed 
case  3.2.  “Frozen”  ions  indicates  particles  were  placed  in  the  problem  by  the  code  in  a  cold  beam.  The  right 
electrode  size  is  in  multiples  of  the  left  electrode. 
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Figure  3:  Vx-x  phase  space  plots  for  lOOOeV  “filling”  beam.  Left:  Current  neutralized,  case  1.3.  Right:  Density 
neutralized,  case  1.4.  Top:  Ion  phase  space.  Bottom:  Electron  phase  space.  t=1.2e-6  s 


REFERENCES 


‘  Ramo-W ooldridge  Staff.  “Electrostatic  Propulsion.”  Proceedings  of  the  IRE,  Vol.  48,  No.  4,  April  1960. 

French,  Park.  “Circular  Beam  Neutralization.”  Progress  in  Astronautics  and  Rocketry  Volume  5:  Electrostatic 
Propulsion,  p.  237-250,  Academic  Press,  New  York,  1961 . 

3  Mirels,  H.  “On  Ion  Rocket  Neutralization.”  Progress  in  Astronautics  and  Rocketry  Volume  5:  Electrostatic 
Propulsion,  p.373-381,  Academic  Press,  New  York,  1961. 

Seitz,  R.N.  et  al.  Present  Status  of  the  Beam  Neutralization  Problem.”  Progress  in  Astronautics  and  Rocketry 
Volume  5:  Electrostatic  Propulsion,  p.383-422,  Academic  Press,  New  York,  1961. 

Buneman,  O.  and  Kooyers,  G.  Computer  Simulation  of  the  Electron  Mixing  Mechanism  in  Ion  Propulsion.” 
AIAA  Journal,  Vol.  1,  No.  1 1,  p.2525-2528,  November  1963. 

6  Wadhwa,  R.P.  et  al.  “Two-Dimensional  Computer  Experiments  on  Ion-Beam  Neutralization.”  AIAA  Journal  Vol. 
3,  No.  6,  p.1076-1081,  June  1965. 

Tajmar,  M.  and  Wang,  J.  “Three-Dimensional  Numerical  Simulation  of  Field-Emission-Electric-Propulsion 
Neutralization.”  Journal  of  Propulsion  and  Power,  Vol.  16,  No.  3,  p.536-544,  May  2000. 

Tajmar,  M.  and  Wang,  J.  “Three-Dimensional  Numerical  Simulation  of  Field-Emission-Electric-Propulsion 
Backflow  Contamination.”  Journal  of  Spacecraft  and  Rockets,  Vol.  38,  No.  1,  p.69-78,  January  2001. 

9Pawlik,  E.V.  “Neutralization  of  a  Movable  Ion  Thruster  Exhaust  Beam.”  JPL  Space  Programs  Summary  37-58 
Vol.  Ill,  1969. 

10  Wang,  et  al.  “Three-Dimensional  Particle  Simulations  of  Ion  Propulsion  Plasma  Environment  Deep  Space  1.” 
Journal  of  Spacecraft  and  Rockets,  Vol.  38,  No.  3,  May  2001 . 

11  Othmer,  C.  et  al.  “Three-dimensional  simulations  of  ion  thruster  beam  neutralization.”  Physics  of  Plasmas  Vol 
7,  No.  12,  p.5242-5251,  Dec.  2000. 

Othmer,  C.  et  al.  “Numerical  Simulation  of  Ion  Thruster-Induced  Plasma  Dynamics  -  The  Model  and  Initial 
Results.”  Advances  in  Space  Research,  Vol.  29,  No.  9,  p.1357-1362,  Elsevier  Science  Ltd.,  UK,  2002. 

Othmer,  C.  et  al.  “Numerical  Parameter  Studies  of  Ion-Thruster-Beam  Neutralization”  Journal  of  Propulsion  and 
Power,  Vol  19,  No.  5,  p.953-963,  September-October  2003. 

14  Verboncoeur,  J.P.,  M.V.  Alves,  V.  Vahedi,  and  C.K.  Birdsall,  “Simultaneous  Potential  and  Circuit  Solution  for  Id 
bounded  Plasma  Particle  Simulation  Codes,”  J.  Comp.  Physics,  104,  pp.  321-328,  February  1993. 

15  Vahedi,  V.,  C.K.  Birdsall,  M.A.  Lieberman,  G.  DiPeso,  and  T.D.  Rognlien,  “Verification  of  frequency  scaling 
laws  for  capacitive  radio-frequency  discharges  using  two-dimensional  simulations,”  Phys.  Fluids  B  5  (7)  pp  2719- 
2729,  July  1993. 

16  V.  Vahedi  and  G.  DiPeso,  “Simulataneous  Potential  and  Circuit  Solution  for  Two-Dimensional  Bounded  Plasma 
^Simulation  Codes,”  J.  Comp.  Phys.  131,  pp.  149-163, 1997.  (Work  begun  at  U.C.  Berkeley) 

17  Gatsonis,  N.A.  and  Spirkin,  A.  “Unstructured  3D  PIC  Simulations  of  Field  Emission  Array  Cathodes  for 
Micropropulsion  Applications.”  38th  Joint  Propulsion  Conference,  Indianapolis,  IN,  July  7-10  2002. 

18  Spirkin,  A.  and  Gatsonis,  N.A.  Unstructured  3D  PIC  Simulation  of  Plasma  Flow  in  a  Segmented  MicroChannel. 
36th  AIAA  Thermophysics  Conference,  Orlando,  FL,  June  23-26  2003. 


